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0\ ■ Abstract 

Self-consistent system of spinor, scalar and BI gravitational fields is considered. Exact solutions 
to the field equations in terms of volume scale of the BI metric are obtained. Einstein field 
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' equations in account of the cosmological constant A and perfect fluid are studied. Oscillatory 
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■ mode of expansion of the universe is obtained. It is shown that for the interaction term being a 
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. power function of the invariants of bilinear spinor forms and A > and given other parameters, 
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e.g., coupling constant, spinor mass etc., there exists a finite range of integration constant which 
generates oscillatory mode of evolution. 
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1. INTRODUCTION 

The discovery of the cosmic microwave radiation has stimulated a growing interest 
in anisotropic, general-relativistic cosmological models of the universe. The choice of 
anisotropic cosmological models in the system of Einstein field equation enable us to 
study the early day universe, which had an anisotropic phase that approaches an isotropic 
one BI universe is the simplest model of an anisotropic universe and eventually 

evolves into a Freidmann- Robertson- Walker (FRW) universe Q, if filled with matter obey- 
ing p = C^e, < 1, where e and p are the energy and pressure of the material field, 
respectively. Since the present-day universe is surprisingly isotropic, this feature of the BI 
Universe makes it a prime candidate for studying the possible effects of an anisotropy in the 
early universe on present-day observations. 

In this paper we study the self-consistent system of spinor, scalar and BI gravitational 
fields in presence of perfect fluid. Solutions of Einstein equations coupled to a spinor and 
a scalar fields in BI spaces have been extensively studied by Saha and Shikin 0, 0, 0, ■ 
In those papers, the field equations were solved qualitatively. In this report, we consider 
some key equations occurred in those papers. Initial value problem has been posed for those 
equations and solved numerically. Application of numerical methods enables us to view this 
problem from totally different angle and gives rise to some interesting results previously 
unknown. 

2. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 

The action of the nonlinear spinor, scalar and gravitational fields can be written as 



where R is the Ricci scalar and L is the spinor and scalar field Lagrangian density chosen 
in the formlsf 



Here A is the coupling constant and F is some arbitrary functions of invariants generated from 
the real bilinear forms of a spinor field. We choose F to be the function of / = S*^ = {ipipY 




(2.1) 
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(2.2) 
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and J = = {i'4'_T''4'Y -I i-e., F = F{I,J), that describes the nonhnearity in the most 
general of its form p. As one sees, for A = we have the system with minimal coupling. 
The gravitational field in our case is given by a Bianchi type I (BI) metric in the form 

ds'^ = alidx^f - aj^dx^f - al{dx^f - al{dx^f, (2.3) 

with oq = 1, x*^ = ct and c = 1. The metric functions {i = 1,2,3) are the functions of 
time t only. 

Variation of 1)2.11) with respect to spinor field ip {ip) gives nonlinear spinor field equations 



iYV^^tfj -mtfj + Vtlj + gi-f^^ = 0, (2.4a) 
iV ^ipY + rmp - Vi> - gi4>-f^ = 0, (2.4b) 

where we denote 

V = XSip^^^'^'dF/dl, g = XPip^^^'^'dF/dJ, 
whereas, variation of ()2.1)1 with respect to scalar field yields the following scalar field equation 



d 



(v^^7^^(l + AF)^,^) =0. (2.5) 



■s/—g dx^ 

Varying 1^2.1^ with respect to metric tensor g^iy one finds the gravitational field equation 
which in account of cosmological constant A has the form 

(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 

Here k is the Einstein gravitational constant and over-dot means differentiation with respect 
to t. The energy-momentum tensor of the material field is given by 
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Here = {£, —p, —p, —p) is the energy-momentum tensor of a perfect fluid. Energy 

e is related to the pressure p by the equation of state p = C,e. Here C, varies between the 
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interval < ( < 1, whereas C = describes the dust Universe, C = I presents radiation 
Universe, | < C < 1 ascribes hard Universe and ( = 1 corresponds to the stiff matter. 
The Dirac matrices 'y^{x) of curve space-time are connected with those of Mincowski space 
as 

7^ = 77a^, 7m = 7«m. = 0,1, 2, 3. (2.8) 
The explicit form of the covariant derivative of spinor is j3| 

V^2P = d^i:-T^i:, V^iJ = d^ij + i:T^, /i = 0,l,2,3, (2.9) 

where r^(x) are spinor affine connection matrices. For the metric ()2.3|1 one has the following 
components of the affine spinor connections 

= (l/2)d^7'^7°- (2.10) 

We study the space-independent solutions to the spinor and scalar field Eqns. ()2.4|) and 
()2.5|) so that ijj = tpit) and ip = ip{t). defining 

r = aQaia2a3 = \J —g (2-11) 

from ()2.5|1 for the scalar field we have 

^ = C [ [r{l + \F)]-^dt. (2.12) 



Setting Vj{t) = ^^Ijjit), j = 1,2,3,4, in view of ((231) and ^TW) from f^:^ one 
deduces the following system of equations: 

Vi + i{m ~ V)Vi - gVs = 0, (2.13a) 

V2 + iim - V)V2 - = 0, (2.13b) 

V3 - iim - V)V3 + QVi = 0, (2.13c) 

V4 - iim - V)V4 + gV2 = 0. (2.13d) 

/^From ()2.4a|) we also write the equations for the bilinear spinor forms S, P and = 

^o-2^A° = 0, (2.14a) 

Po -2(771-1)) A° = 0, (2.14b) 

i° + 2(777 -r')Po + 2^5o = 0, (2.14c) 
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where Qo = tQ, leading to the relation S'^ + P"^ + {A^f = C'^/t'^, = const. As 

one sees, for F = F{I) ()2.14a|) gives 5* = Cq/t, whereas for the massless spinor field with 
F = F{J) ()2.14bjl yields P = Dq/t. In view of it for F = F{I) we obtain the following 
expression for the components of spinor field 

(2.15) 

with Ci being the integration constants and are related to Cq as Cq = Cf + C2 — ~ Cf. 
Here (3 = Jim- V)dt. In case of F = F{J) for the massless spinor field we get 

(2.16) 

The integration constants Di are connected to Dq hj Dq = 2 (Df + D^ — D^ — D|). Here we 
set cr = / Qdt. 

Once the spinor functions are known explicitly, one can write the components of spinor 
current = ^'y'^if)^ the charge density of spinor field q = (jo-j^Y^'^, the total charge of spinor 



field Q = J g^/ —^gdxdydz, the components of spin tensor S'^'^''' = ^ij) [^^"^ a^" + a'^'^Y}'^ 
—00 

and other physical quantities. 

Let us now solve the Einstein equations. In doing so we first write the expressions for 
the components of the energy-momentum tensor explicitly: 



TO = mS + Cy2T^{l + XF)+e, 

= T| = Ti = VS + gP- C72r2(l + XF) - p. 
In account of ()2.17j) from ()2.(iaj) . ()2.6bj) . ()2.(icj) we find the metric functions 



(2.17) 



X, / [Tit)]-'dt 



ai{t) = Ar'/^exp 

with the constants of integration Di and Xj obeying 

3 3 

1=1 1=1 



1,2,3, (2.18) 
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As one sees from fj2.18|) for r ~ with n > 1 the exponent tends to unity at large 
t, and the anisotropic model becomes isotropic one. Let us also write the invariants of 
gravitational field. They are the Ricci scalar Ii = R ^ l/r, I2 = R^uR^^ = R^^Ry ~ 
and the Kretschmann scalar Is = Rap^^R'^^^'^ ^ As we see, the space-time becomes 

singular at a point where r = 0, as well as the scalar and spinor fields. Thus we see, all the 
functions in question are expressed via r. In what follows, we write the equation for r and 
study it in details. 

Summation of Einstein Eqns. ()2.6a|) . ()2.6b|) . ()2.6c|) and ()2.6d|) multiplied by 3 gives 



r _ 3 
r ~ 2 



K (mS + VS + gP + e -p^ - 3A. (2.19) 



From energy-momentum conservation law T^.^ = 0, in account of the equation of state 
p = (e, we obtain 

e = eo/T'+<, p = Ceo/T'+<. (2.20) 

In our consideration of F as a function of /, J ot I ± J we get these arguments, as well as 
V and Q as functions of r. Hence the right-hand-side of ()2.19p is a function of r only. In 
what follows we consider the case with F = F{I). Recalling the definition of T>, in view of 
fjTT^ and ^nn^ the Eqn. (jTTT?|) can be written as 

f = ^(r,p), (2.21) 

where we denote 

^ ^^(^rnCo + XC',C'Fj{t)/t'{1 + AF(r))2 + eo{l - 0/t<) - 3Ar , (2.22) 

and p = {k, X,m,CQ,C,eQ,(,-A} is the set of the parameters. Here we take into account 
that S = Co/t. From mechanical point of view the Eqn. ()2.21|) can be interpreted as an 
equation of motion of a single particle with unit mass under the force JF(r, p). Then the 
following first integral exists 



f = ^/E-UiT,p). (2.23) 
Here E is the integration constant and 

K^mCoT - CV2(1 + XF) + eoT-'^^ - Ar^ 
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is the potential of the force T . We note that the radical expression must be non-negative. 
The zeroes of this expression, which depend on all the problem parameters p define the 
boundaries of the possible rates of changes of rii). Note that setting m = in ()2.22j] we 
come to the case of massless spinor field with F = F{J) or F = F{I ± J). 

We formulate the initial value problem for the Eqn. 1)2.211) with initial condition 



which we solve numerically. 

Since r is the volume-scale, it cannot be negative for every t > 0. On the other-hand, BI 
space-time models a non-static universe, i.e., the derivative f{t) should be nontrivial at the 
initial moment t = 0. This leads to to fact that for fixed p, the constant E and the initial 
value of T are inter-related in the sense that for a given E the value of Tq should belong to 
some interval. 

In what follows, we numerically solve ()2.21)1 for some concrete form of F. 
Let us choose F as a power function of S*, namely, F = S". In this case setting Cq = 1 
and C = 1 we obtain 



Note that the nonnegativity of the radical in ()2.23|) in view of f|2.25|) imposes restriction on 
r from above in case of A > 0. It means that in case of A > the value of r runs between 
and some Tmax, where Tmax is the maximum value of r for the given value of p. This equation 
has been studied for different values of parameters p. Here we demonstrate the evolution of 
r for different choice of tq for fixed "energy" E and vise versa. 

As the first example we consider massive spinor field with m = 1. Other parameters are 
chosen in the following way: coupling constant A = 0.1, power of nonlinearity n = 4, and 
cosmological constant A = 1/3. We also choose ( = 0.5 describing a hard universe. 

In Fig. ^we plot corresponding potential W(r) multiplied by the factor 2/3. As is seen 
from Fig. ^ and Fig. El choosing the integration constant E we may obtain two different 
types of solutions. For E > 0.5 solutions are non-periodic, whereas for -Emm < E < 0.5 the 
evolution of the universe is oscillatory. 



r(0) = ro>0. 




(2.24) 



with the potential 




(2.25) 
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FIG. 1: Perspective view of the potential W(t). 




FIG. 2: Perspective view of r for different choice of E. 



As a second example we consider the massless spinor field. Other parameters of the 
problem are left unaltered with the exception of (. Here we choose C = 1 describing stiff 
matter. It should be noted that this particular choice of ( gives rise to a local maximum. 
This results in two types of solutions for a single choice of E. 

As one sees from Fig. |21 for E > M there exists only non-periodic solutions, whereas, for 
i?min < E < —0.5 the solutions are always oscillatory. For E G (—0.5, M) there exits two 
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FIG. 3: Perspective view of the potential U{t) with BI universe being filled with stiff matter. 
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FIG. 4: Perspective view of r for different choice of tq with E G (—0.5, M). 

types of solutions depending on the choice of tq. In Fig. E] we plot the evolution of r for 
E G (—0.5, M). As is seen, for tq G (0, A) we have periodical solution, but due to the fact 
that r is non-negative, the physical solutions happen to be semi-periodic. For tq G {B, C) 
we again have oscillatory mode of the evolution of r. This two region is separated by a 
no-solution zone (A, B). 



Case:m = 0, X = 0.1,n = 4, C = 1, A= 1/3 
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3. CONCLUSIONS 

A self-consistent system of spinor, scalar and gravitation fields has been studied in pres- 
ence of perfect fluid and cosmological term A. Oscillatory mode of evolution of the universe 
is obtained. It is shown that for the interaction term being a power function of the invariants 
of bilinear spinor forms, the oscillatory solution is possible if A-term is positive. It is also 
shown that only for a finite range the integration constant E there exists oscillatory mode 
of evolution. It should be emphasized that a third type of solution is possible, if the BI 
universe is filled with stiff matter. 
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